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ABSTRACT 

We explore whether global observed properties, specifically half-light radii, mean surface brightness, 
and integrated stellar kinematics, suffice to unambiguously differentiate galaxies from star clusters, 
which presumably formed differently and lack dark matter halos. We find that star clusters lie on 
the galaxy scaling relationship referred to as the Fundamental Manifold (FM), on the extension of 
a sequence of compact galaxies, and so conclude that there is no simple way to differentiate star 
clusters from ultra-compact galaxies. By extending the validity of the FM over a larger range of 
parameter space and a wider set of objects, we demonstrate that the physics that constrains the 
resulting baryon and dark matter distributions in stellar systems is more general than previously 
appreciated. The generality of the FM implies I) that the stellar spatial distribution and kinematics 
of one type of stellar system do not arise solely from a process particular to that set of systems, such 
as violent relaxation for elliptical galaxies, but are instead the result of an interplay of all processes 
responsible for the generic settling of baryons in gravitational potential wells, 2) that the physics of 
how baryons settle is independent of whether the system is embedded within a dark matter halo, and 
3) that peculiar initial conditions at formation or stochastic events during evolution do not ultimately 
disturb the overall regularity of baryonic settling. We also utilize the relatively simple nature of star 
clusters to relate deviations from the FM to the age of the stellar population and find that stellar 
population models systematically and significantly over predict the mass-to-light ratios of old, metal- 
rich clusters. We present an empirical calibration of stellar population mass-to-light ratios with age 
and color. Finally, we use the FM to estimate velocity dispersions for the low-surface brightness, outer 
halo clusters that lack such measurements. 

Subject headings: galaxies: formation, evolution, star clusters, structure, stellar content 



1. INTRODUCTION 

Stellar systems are conventionally divided into two 
broad categories, stellar clusters and galaxies, depending 
on whether they lie within a dark matter halo. Yet this 
scheme faces conceptual and practical difficulties. For ex- 
ample, the dynamical masses of cert ain "galaxies" (ultra- 
compact dwarf galaxies or UCDs; Mieskc ct al.l 12008ft 
suggest little or no dark matter, as does the modeling of 
tidal dwarfs ([Barnes fc Hcrnquistl 119921 ). Furthermore, 
even for indisputable galaxies like normal ellipticals, it 
can be difficult to infer dark matter fro m kinematic mea- 
surem ents within the luminous radius fCappellar i et al.1 
2007). It is the spatial concentration of baryons rel- 
ative to dark matter within the luminous region that 
distinguishes UCDs from other galaxies and that makes 
dark matter easier to measure in spirals than ellipticals. 
Rather than differentiating stellar systems with a binary 
classification — dark matter halo vs. no halo — that is 
difficult to measure, we consider baryon vs. dark matter 
concentration, which is easier to observe and allows for 
a continuum of properties. 

Our pr eviou s work on scaling relations (Zaritsky et al. 
I2006alf b1. 12008| ) shows that galaxies lie on a progression 
of star formation efficiencies and baryon concentrations 
within dark matter halos. The mass-to-light ratio T e 
within the half-light radius r e is one gauge of the rel- 
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ative concentration of baryons to dark matter. The 
large T e 's of dwarf spheroidals (dSphs) reflect the im- 
portance of dark matter within their luminous region, 
while the lower T e 's of ellipticals arise from the dom- 
inance of the baryons there. One could hypothesize a 
more extreme case than that for ellipticals, where the 
baryons are packed so tightly that the dark matter is dy- 
namically negligible within r e . It would be difficult to 
differentiate such a low T e galaxy, if it were low mass, 
from a star cluster, which is devoid of dark matter. I f 
they have any dark matter, UCDs (jMieske et al.l 12008ft 
may be examples of galaxies whose spatial and kinematic 
structure thus resemble those of star clusters. Are there 
other observables that reveal their presumably distinct 
formation and evolutionary paths? Remarkably, compa- 
rably low mass systems can have low T e (stellar clusters, 
UCDs) or high T e (dSphs) . What variations in dark halo 
properties of baryonic processes explain this wide range 
of baryon to dark matter concentrations? 

Such questions are fundamental to understanding the 
nature of low mass stellar systems, the basic characteris- 
tics of galaxies, and the ways in which baryons can settle 
into a dark matter halo. Here we u se a new scaling rela- 
tion (the Fu ndamental Manifold of iZaritskv et al.f (FM: 
I2006all2008ft ) to explore the first question above, whether 
star clusters can be distinguished from low mass galaxies 
with similar T e 's and the implications for their evolution. 
In certain projections of paramete r space, globular 



cfuslcrs are distinc t from galaxies (|Kormendvl 119851 : 
iBurstein et al.l fl997h . (Then again, in certain projec- 
tions, dwarf galaxies are distinct from giant galaxies.) 
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Projections of distributions in multi-dimensional param- 
eter space can be misleading. It is critical that the de- 
fined parameter space 1) has physically meaningful axes, 
so that one can interpret the results, 2) produces a galaxy 
distribution as thin as that allowed by observational er- 
rors (i.e., that has N — 1 dimensionality), thus ensuring 
that the full dimensionality of the distribution is cap- 
tured, and 3) exploits the full range of galaxy properties 
to maximize the discriminatory power of the observa- 
tional space. In other words, we seek the minimum set of 
physically-meaningful parameters that captures the ma- 
jority of the variance in galaxy properties. 

The Fundamental Manifold, so named in ref- 
erence to its ante c edent , the Fundamental Plane 
piorgovski fc Davis! 119871 ; iDressler et all 11987ft . is 
tightly populated by galaxies of all morphological types 
and luminosities. By asking where star clusters lie rela- 
tive to this manifold, we have a means to quantify how 
distinct star clusters are from galaxies. The relationship 
is straightforward 

log r e = 2 log V - log I e - log T e - C, (1) 

where r e is the half-light radius in kpc, V represents the 
internal velocity of the system and is a combination of 
the velocity dispersion, a, and the rotational velocity, v r , 
in km s^ 1 , I e is the mean surface brightness within r e in 
solar luminosities per square parsec, T e is the total mass- 
to-light ratio of the galaxy within r e in solar units, and 
C is a normalization constant. Additionally, and what 
leads to the finding that the systems populate a 2-D man- 
ifold, is the empirical determination that T e is primarily 
a function of V and I e . The functional form of T e (V, I e ) 
is not simple, hence the broad term manifold, and while 
the data used to derive the fitted surface extend over 
many orders of magnitude along each of the relevant pa- 
rameters, the function is unconstrained over much of the 
plausibly accessible parameter space and must be used 
with caution. For ex ample, the function as described by 
iZaritskv et ail (2008) was unconstrained in the region of 
parameter space occupied by star clusters. 

Equation (1) originates from the virial theorem, and 
when combined with the finding that T e = T e (V,I e ), 
it limits the region of parameter space occupied by the 
galaxy population. In other words, the virial theorem 
alone allows infinitely more possible configurations that 
those found to populate the FM. As such, the FM places 
a strong constraint on galaxy formation and evolution 
models. The simplicity of Equation (1) and its relation 
to the virial theorem satisfies our first condition above. 
The empirical finding that the scatter about this relation 
is small for the most diverse set of galaxies available, 
and yet captures their variety in kinematics, size, and 
luminosity, demonstrates that we have included the full 
dimensionality of the family of galaxies, at least to within 
the observed scatter, and therefore satisfies our second 
and third conditions. 

Before testing whether star clusters lie on the FM, it is 
important to understand what the relationship implies 
about galaxies. First, satisfying Equation (1) requires 
that galaxies be in virial equilibrium and that they have 
what we will refer to as "structural terms" that have 
a small degree of scatter. Hereafter, we use the term 
structure to encompass both the spatial distribution and 
kinematics of stars and dark matter. These structural 



terms are the calculated numerical coefficients of the v 2 
and GM/r terms one would obtain when evaluating the 
kinetic or potential energy terms in the virial theorem. 
For example, an imaginary galaxy of constant density 
out to r e would have a different total potential energy 
than a galaxy that consists of the same amount of mass 
in a ring of radius r e , although the potential energy of 
either galaxy can be expressed as being proportional to 
GM/r e . It is impossible to numerically tally the terms 
in the virial theorem because one never has complete in- 
formation regarding the phase space distributions of the 
stars and dark matter in a galaxy. Nevertheless, the low 
scatter and lack of systematic differences among differ- 
ent galaxy types when applying Equation (1) indicates 
that the structural terms, which we have blithely incor- 
porated into the constant term, are not detectably dif- 
ferent among galaxies. Second, and perhaps somewhat 
more surprising, is that T e is principally determined by 

V and I e , and therefore only modestly related to other 
factors that might have had an influence on T e , beyond 
whatever joint influence those factors may have had on 

V and I e . As such, there is a strong connection between 
the overall structure of a galaxy and the mass-to-light 
ratio within r e . It is this second part in particular, the 
interdependence of V, I e and T e , that one might suspect 
would be qualitatively distinct for star clusters if their 
formation or evolution mechanism(s) is wholly different 
than that of galaxies. Consider that dSphs, with their 
< 10 km s _ velocity dispersions, are quite different in 
appearance from star clusters with similarly low velocity 
dispersion. 

Our plan is to test whether star clusters can be incorpo- 
rated easily into the FM construct. If they can be, then 
they can be described as a low-mass, high-concentration 
extension of the family of galaxies because the lack of 
dark matter within r e has produced no detectable differ- 
ence in their structure relative to galaxies. On the other 
hand, if clusters lie significantly off the FM then their 
structure is sufficiently different from that of galaxies to 
suggest differences in their formation and/or evolution. 
To provide a specific example, one could imagine a sys- 
tem like a star cluster that is instead dominated by dark 
matter (consider turning some fraction of the stars into 
dark matter). This system would have the kinematics 
and size typical of normal star clusters, but a very dif- 
ferent surface brightness and mass-to-light ratio. If we 
then allow the dark matter to follow a different radial 
density profile, all the observables would be affected and 
we would not necessarily expect this system to lie on the 
same relationship between T e , V , and I e . The converse 
does not hold in the sense that systems on the FM do 
not necessarily have the same evolutionary history. 

Although a better understanding of the origin of galax- 
ies and star clusters is the primary goal, universal scal- 
ing relationships provide numerous ancillary benefits. 
For example, we are familiar with the distance deter- 
minations using the Fundamental Plane (FP) and Tully- 
Fisher (TF) relationships. Unlike those scaling relations, 
the FM also provides a measure of the total mass-to-light 
ratio within r e . Given that star clusters lack dark matter 
halos, we will use the FM to constrain the stellar mass- 
to-light ratio, Y», as a function of cluster parameters 
such as main sequence turn-off age and color. We will 
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then compare these measures to values calculated using 
stellar population models as a test of those models. 

In summary, by considering star clusters and UCDs, 
we test whether the FM formalism extends to the low- 
mass, stellar-dominated extremum of gravitationally- 
bound stellar systems in §3. We find that star clusters do 
lie on the FM but with apparently greater scatter than 
observed among the galaxies. We then explore both ob- 
servational and physical sources of that scatter. We find 
that there is a correlation between cluster age and devi- 
ation from the FM, in the sense expected for variations 
in T e due to the evolution of the stellar population, but 
in quantitative disagreement with the predictions from 
stellar population models. We explore this behavior and 
present empirical derivations of the relationships between 
T e , age, and B — V in §4. We find that the dominant 
source of scatter in the current data set is observational, 
arising from the uncertainties in the cluster velocity dis- 
persions. We use the FM to make predictions for the 
velocity dispersions of those clusters without measured 
dispersions. In §5 we summarize our findings and dis- 
cuss the implications. In an Appendix we revisit previ- 
ous claims that are in apparent conflict w ith our results 
(jForbes et al.ll2008t iTollerud et aLll2010f ) and provide a 
resolution. 

2. THE DATA 

We draw the measurements of structural parame- 
ters, velocity dispersions, modeled stellar population 
mass-to-light ratios (T* imo( j), ages, colors (B — V), 
and metallicities ((Fe/H)) from the compilation by 
iMcLaughlin fe van der Marell (|2005f ) of Local Group star 
clusters. That sample consists of 153 spatially resolved 
clusters from a range of nearby galaxies. Within that 
sample there is a subset of 57 clusters for which they 
have also compiled velocity dispersion measurements. As 
such, it includes clusters of all ages, a wide range of 
metallicities, and the full range of environments found 
in the Local Group. 

The data for stellar systems ranging from dwarf 
spheroidals to galaxy cluster spheroids (CSphs to dis- 
tinguish them from stellar clusters th roughout, see 
IZaritskv et all [2001 are drawn from Zaritskv et al.l 
( 2008T) and for ultra-compact dwarfs from Micsk ~t al.l 
(2008), selecting only objects with estimated masses 
> 5 x 10 6 M Q to distinguish the UCDs from the stellar 
clusters in their sample. 

3. STAR CLUSTERS AND THE MANIFOLD 

To define the manifold, we use only t he sp heroidal 
galaxies and CSphs from IZaritskv et al.l (|2008h as our 
comparison sample to minimize potential problems aris- 
ing from composite stellar populations and to limit our- 
selves to systems supported by velocity dispersion rather 
than rotation. In the vast majority of local galaxies, ei- 
ther v r or a dominates, and therefore we have not yet 
truly tested whether V 2 = 0.5v 2 + a 2 is preferable to 
other parameterizations, such as V 2 = max[0.5v 2 , a 2 ], 
or whether 0.5 is the optimal coefficient for the v r 
term (this topic will be addressed in another paper, 
but values of 0.3 to .5 are plausible in idealized mod- 
els; iWeiner et al.ll2006h and such parameter i zations have 
been used elsewhere (iBurstein et al.l 119971 ; iKassin et all 
2007; Covington et al. 2010). Here, we avoid questions 



regarding the incorporation of rotation supported sys- 
tems in the FM formalism by only analyzing systems 
in which a dominates. The low relative rotation of 
Milky Way clusters is inf erred from their low ellipticitics 
(jFreeman fc Norrisl [l98l . This argument is weaker for 
other cluster populations, such as those in the Magellanic 
Clouds, where larger ellipticities are sometimes found. 
However, even in the C louds, the cluster ellip ticities are 
generally less than 0.5 (jHill fc Za ritskv 2006), similar to 
the ellipticity of one of the most extreme MW clusters, 
to Cen, which has v/a ~ 0.5 (jReiins et al.ll2006lh There- 
fore, for the remainder of this study, we adopt V = a. 
Using only spheroidal galaxies further simplifies our work 
by removing systems with ongoing or recent star forma- 
tion, for which T e may have a dramatic dependence on 
mean stellar age. 

The final relevant technical issue with the galaxy sam- 
ple is the heterogeneous nature of the source samples 
and how one should compare it to the cluster sample. 
In particular, the cluster sample contains U-band pho- 
tometry, while the galaxy sample contains a range of 
photometry from V to /-bands. Equation (1) is band- 
independent in that I e and T e for a particular object 
need to be in the same band, but that band can be 
different for different objects. However, the manifold 
that we fit to describe T e , T e jn, is band-specific and 
has been derived using the heterogeneous galaxy sam- 
ple, corrected to the /-band. The values of I e used to 
determine T e jn are in solar units, so the filter depen- 
dence arises from the difference in colors between the 
galaxy and the Sun when converting L e to I e . In ef- 
fect we are assuming solar colors for the galaxies, al- 
though any uniform color difference among all galaxies 
is absorbed into the calibration of C in Equation (1). 
The remaining concerns are color differences that vary 
across the sample and band- related differences in r e . If 
these exist, then they are inappropriately incorporated 
into our fit for T e . One clear path to progress is obtain- 
ing a single, uniform data set for the further study of the 
FM. However, in the interim, we proceed with the mot- 
ley sample available and caution that for systems with 
colors much different than solar colors there will be a 
yet unknown color-correction to T e ji t . The empirically- 
confirmed small scatter amo ng galaxies, even wh en in- 
cluding star forming galaxies (jZaritskv et al.| [2008) . sug- 
gests that this effect will be negligible for a sample con- 
sisting only of spheroidal galaxies, but could be a factor 
among the star clusters, which span a wider age range. 

We revaluate T e jj t using only th ose galaxies, the 
spheroids, from IZaritskv et al.l (|2008l ) that are a dom- 
inated, including galaxy cluster spheroids (CSph). In 
doing so, we have also reevaluated the const ant C in 
Equation (1) using the ICappellari et al.l ([2007D sample, 
for which they have full dynamical modeling fit to 2- 
dimensional spectral data for spheroidal galaxies, and 
find C = 0.75. In Figure Q] we show where the 57 glob- 
ular clusters with measured velo city dispersions and the 
UCDs with M > 5 x 1O 6 M from lMieske "eTall (pOOl ) fall 
when using this fitting function, which was not derived 
using either the star clusters or UCDs. In addition, we 
show with small points the range of possible locations 
in this projection of the extrapolated relationship cov- 
ered by random, < lcr deviations of each of the fitted 
coefficients. The star clusters and UCDs lie within the 
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Fig. 1. — The edge-on view of the FM for spheroidal stellar systems, including star clusters and ultra-compact dwarf galaxies (UCDs). 
We are extrapolating the fit from the region populated by systems with logr e > —1 (i.e. excluding both UCDs and stellar clusters). We 
illustrate the statistical uncertainties in that extrapo lation with dots representing projections of the surface where individual fit coefficients 
are altered by lc. Spheroidal stellar systems from Zaritsky et al. (2008) arc plotted as filled black circles, clusters with tabulated a's 
HMcLaughlin fc~a n dor Marel 2005, and references therein) are plotted as filled red circles, and UCDs with M > 5 X 10 5 Mq from 
IMieske et £01 12008') are plotted as filled blue squares. The extrapolated fit, within its lcr uncertainties, passes through the majority of 
the stellar clusters and UCDs, although the scatter among the clusters is g reater than among the galaxies. Uncertainties are not plotted, 
but the dominant uncertai nty in the cluster pop ulation is discussed in i|3.2l and the uncertainties in the comparison galaxy population are 
comparable to the scatter (Zaritsky et al. 2008). The clusters that deviate well below the fit are the youngest clusters. 



extrapolated version of T e ,fit given the fitting uncertain- 
ties. Because the extrapolated fitting function works well 
in the mean for stellar clusters, we infer that clusters are 
indeed the low-mass, limiting case of high-stellar-mass- 
fraction systems and refit T e jit using the cluster and 
UCD data (Figure d]). The resulting fit is 



T eJlt = 1.49 - 0.32 log V - 0.83 log I e 
+ 0.24 log 2 ^ + 0.12 log 2 I e 
-0.02 log V log I e 



(2) 



and is used hereafter. 

Projections of this function along the various axes (Fig- 
ure [3]) illustrate how difficult it is to see the multidimen- 
sional nature of the distribution and how certain popula- 
tions of objects appear to stand out as different despite 
populating the same surface. Because there is no phys- 
ical motivation for this particular form, the addition of 
new data, particularly those that probe a new part of 
parameter space, could require changing the functional 
form. However, the implications of the FM are indepen- 



dent of the particular functional form, depending only 
on the finding that T e ji t is principally a function of V 
and I e . The apportionment of the constant term between 
T e ji t and the C in Equation (1) is set by the external 
calibration to independently measured values of T e and 
by construction in Equation (1) that sets the two sides to 
be equal. The scatter about the FM is evidently larger 
for the star clusters than it is for the galaxies (0.35 vs. ~ 
0.1; lZaritskv et al.| [2008). so we turn to discussing possi- 
ble sources of that scatter. 

There are various potential causes of the larger scat- 
ter, including structural variations among clusters that 
are not present in galaxies, proportionally larger uncer- 
tainties in at least one of the observed parameters for 
the clusters, and the extrapolation of a fitted function 
beyond the parameter range over which it is constrained. 
Among the observed parameters, the likely source of 
scatter is V, which is more poorly measured in these 
low-velocity systems than in normal galaxies, and so we 
will begin our discussion by assessing the impact of the 
quoted uncertainties in V. However, of greater interest 
is the possibility that some of the scatter is driven by a 
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Fig. 2. — The edge-on view of the FM for spheroidal stellar sys- 
tems, where the empirical fit for T e f it includes both the stellar 
clusters and UCDs. Symbols are as described for Figure [l] 

parameter that has not been included in our expression 
for T e jit- This type of effect is often refereed to as a 
"second parameter" question in correlations, although in 
our case it would be a third parameter. If there is such 
a parameter, it has been inadvertently neglected in our 
formulation of the fitting function for T e jn. As such, T e 
as evaluated using Equation (1), T e< FM, should correlate 
with this yet unidentified parameter (as it does with V 
and I e ). 

In closing, we note that there a re two studies 
(jForbes et all 120081; iTollerud et aLll2010H that appear to 
conflict with our claim that stellar clusters satisfy galaxy 
scaling relations. We describe our interpretation of those 
studies in Appendix A and conclude that neither study 
is truly in conflict with the results presented here. 

3.1. Observational Source of Apparent Scatter: Velocity 
Dispersions 

For these low velocity dispersion systems, the uncer- 
tainties in a are fractionally the largest of all the input 
parameters. We now explore the level of scatter about 
the FM expected simply due to the quoted uncertain- 
ties. We use the 51 old (log(age) > 10 Gyr) clusters with 
published velocity dispersion measurements and associ- 
ated uncertainties to avoid any possible effect of age on 
the scatter. We place each cluster exactly on the FM 
by adopting a value for T e that satisfies Equation (1). 
We then draw a thousand cluster samples using Gaus- 
sian distributed errors on the velocity dispersion, using 
the corresponding value for the uncertainty in a for each 
particular cluster. In Figure [4] we plot the distribution 
of the simulated scatter for the thousand trials and the 
observed scatter for comparison. The observed scatter is 
somewhat larger than the typical simulated one, but we 
cannot exclude the possibility that the observed scatter 
is directly due to the a uncertainties with greater than 
80% confidence. We conclude that the observed scatter 
in Figured] is principally due to the uncertainties in V, 
but cannot conclude that no physical sources of scatter 
remain. We test for physical sources of scatter in the 
next three subsections. 

3.2. Physical Sources of Apparent Scatter: Age 

There is at least one physical cause of scatter that 
must play a role, that we can quantify easily, and that 



Fig. 3. — Projection of T e fit along the three observable param- 
eters. The plume of objects rising in the left of the first and third 
panels are the dSphs. The objects near the bottom left of the first 
and third panels are the star clusters and UCDs. It is difficult in 
projection to determine whether a single surface fits all the data 
and, if not, which population to fit. We find that a single surface 
does fit all the objects. 
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Fig. 4. — Distribution of simulated FM scatter for star clusters. 
We assign each cluster a value of a that places it exactly on the 
FM. We then use the quoted uncertainty in a to create random 
samples. The plot shows the resulting scatter distribution for 1000 
samples. The dashed line shows the measured scatter. Although 
the observational scatter is larger than the typical simulated one, 
we cannot confidently exclude the possibility that it is only the 
result of the fractionally larger uncertainties in noBS f° r these low 
a systems. 

we can model — variations in the stellar population 
mass-to- light ratio, T*. In Figure [5] we demonstrate 
that the observed deviations from the FM for star clus- 
ters are, as expected, different for the young and old 
stellar systems. We fit a linear relation in log (age) to 
describe the residual (0.272(log(age) - 10.11) + 0.044) 
and subsequently use this fit and th e cluster ages from 
iMcLaughlin fc van der Marell (|2005l ). obtained from re- 
solved stellar main sequences, to correct T e ,/it for this 
systematic age deviation. Importantly, the residual is 
small for the older clusters, for which the calibration to 
galaxies is expected to be most appropriate. There re- 
main, however, at least two open questions. First, the 
scatter remains large for both young or old clusters (al- 
though, at least for the old clusters, we demonstrated 
above that observational uncertainties can account for 
the scatter). Second, the mean trend only testifies to 
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some age-related effect, not that this effect must be re- 
lated to T«. There are three reasons that the trend is 
likely due to the evolving stellar population rather than 
to age-dependent dynamical effects: 1) the stellar popu- 
lation must, to some degree, evolve and affect T e because 
T e is related to T„, 2) the fitting function works best (i.e. 
the residual is ~ zero) for the oldest clusters, for which 
the calibration for spheroidal galaxies and their similar 
stellar populations should be most appropriate, and 3) 
dynamical effects, such as those involving two-body re- 
laxation, should trend in the opposite direction with age, 
becoming more, rather than less, important in the older 
clusters. Unfortunately, further discussion is limited by 
the lack of clusters with measured a across our entire 
age range, 8 < log(age[yrsl) < 10. Although the entire 
McLaughlin & van der Marcl ( 2005) cluster sample does 
contain clusters in this age range, we need an estimate 
of cr to include them in our discussion. We next describe 
how we use the FM to obtain such estimates. 

Because the right hand side of Equation (1) can be 
written as only a function of V and I e after adopting 
^ejiti we use the tabulated values of I e and r e for the 
clusters to solve for V, which we equate with a (see §3). 
The drawback of this approach is that we have now as- 
sumed, on the basis of Figure [TJ that both Equation (1) 
and Y e jit are applicable to clusters. We account for the 
trend of T e jn with age shown in Figure [5] using the lin- 
ear fit, and then we solve the equation for V. We test this 
procedure on the subset of 57 clusters with measured a in 
Figure El The agreement is encouraging and we proceed 
to estimate a for all 153 clusters. The agreement within 
the uncertainties of the calculated and observed cr's in 
Figure [6] demonstrates that much of the scatter among 
clusters on the FM must come from the uncertainties in 
their respective cr's, as we had surmised on other grounds. 
In other words, we have not created a situation where we 
absorb errors in other quantities or intrinsic scatter into 
our derived values of a. In all subsequent plots, we use 
the inferred it's, but highlight those clusters with mea- 
sured cr's using large open circles to illustrate the degree 
to which any trends are driven by those clusters with 
unconfirmed velocity dispersion estimates. 

An interesting subset of clusters are those of lower 
surface brightness in the outskirts of the Milky Way. 
Because of their distance and naturally low cr, they 
tend to not have a measurements of high relative pre- 
cision. Outside of the compilation of cl uster proper- 
ties bv lMcLaughlin fc van der Marell (|2005t l that we have 
used exclusively, the re are measu rements of a for Pal 5 
(jDehnen et al.ll200l and Pal 14 (jJordi et all [20091. Us- 
ing these values, we find that even these clusters fall 
on the FM relation. As a target for future observers, 
we present the calculated values of cr necessary to place 
similarly unusual clusters on the FM (Table [lj . These 
low dispersions illustrate why high precision spectroscopy 
and a detailed know ledge of the binary populations (cf. 
iDehnen et all (|200l ) are needed to confirm the calcu- 
lated T e 's and to resolve the question of whether such 
systems lie on the FM. 

We now return to investigate the inferred drift in T e vs. 
age by including all star clusters with or without mea- 
sured cr's. We plot the total mass-to-light ratio within r e 
evaluated using Equation (1), T s> fm, vs. age (Figure[7j. 
For those clusters where we estimate cr, the values of T e 
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Fig. 5. — Deviations from the F M vs. age for the clusters in 
McLaughlin & van der Marcl (2005) sample. The dashed line rep- 
resents the best fit line to the relation. As expected, the mean 
deviation is nearly zero for older systems, whose populations are 
most like those of early-type galaxies, from which we derived the 
functional form of T e ^ it . The scatter is larger separately for either 
the old and young populations than for the galaxies, but that is 
mostly, if not wh olly, due to the clusters' proportionally larger cr 
uncertainties (see !|3.H . 
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Fig. 6. — Comparison of a estimated from placing clusters on the 
FM (ofm) to their measured velocity dispersions (o"oss). The 
solid line represents the 1:1 line, while the dashed line represents 
the b est fit line. The procedure works wel l, so we apply it to the 
entire McLaughlin & van der Marel (2005) sample to estimate a 
for those clusters that do not yet have cjobs- 

are not independent of the measured age, and so trends 
should be viewed with caution and compared to those 
obtained from clusters with measured a. Because stellar 
clusters are dark matter free, T e = T». We confirm that 
the value of T e necessary to place clusters on the FM is a 
simple function of age and, as seen in Figure [71 that the 
sense of the trend (higher values of T e for older ages) is 
as expected and matches the behavior determined from 
only the clusters with measured a. Because all but the 
lowest four of the clusters in Figure [Hare old, the appar- 
ent scatter is not the result of age differences within the 
sample. 

3.3. Physical Sources of Apparent Scatter: Structural 

Terms 

Another potential source of scatter is variance related 
to the cluster structural properties. Star clusters can 
evolve dynamically due to their relatively short dynam- 
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Fig. 7. — Total mass-to-light ratio within r e as evaluated using 
Equation (1), T e fm> vs - main sequence turn-off age compiled by 
McLaughlin & van dcr Marcl (2005) using stellar color-magnitude 
diagrams from a set of references. All clusters are plotted using 
the estimated cr's, up m (small open circles). Those clusters that 
also have a measured <r, a"OBS> are highlighted (large open circle). 
For comparison, we recalculate T e} FM using oobs when avail- 
able (filled circles that are offset in age by 0.2 for clarity). On the 
right, we show the distribution of T e fm f° r the old (log(age) >10) 
clusters. The dashed line represents the distribution using apM, 
while the solid line represents the distribution using o-qbs (the 
histograms are scaled arbitrarily for better comparison). The in- 
clusion of all clusters supports the conclusion from Figure \E\ that 
there is a correlation of T e with age. 

ical times (in comparison to the dynamical times for 
galaxies) and therefore should have larger variations in 
their structural terms, which reflect both the stellar kine- 
matics and spatial distribution. We examine whether 
structural differences among the clusters have been in- 
correctly excluded from Equation (f ), or T e ,/it, by look- 
ing for correlations between T e} FM and measurements of 
the structural state of the cluster. In a similar approach, 
McLaughlin (2000) found some evidence for correlations 
between the structural state of Galactic clusters and their 
location relative to the FP. Figure [5] shows the relation- 
ships between T C| j7M and central potential depth, con- 
centration, and galactocentric distance. Whatever trends 
there may be in the lower panels, they are clearly of much 
lower magnitude than that with age, which is evident as 
the bimodal distribution in the upper panels of the Fig- 
ure. We conclude that none of the measures of cluster 
structure that are currently available suggest a relation- 
ship between deviations from the FM and the dynamical 
state of the cluster beyond any directly related to age. 
Because dynamical state depends on more than just age, 
we conclude that stellar population differences, which do 
depend only on age, are the primary driver of the trends 
we observe. 

3.4. Physical Sources of Apparent Scatter: Metallicity 

Next, we search for a dependence of T Bj fm 011 cluster 
metallicity (Figure [9]). There is no evident dependence of 
Y e ,FM on metallicity below {Fe/H} < —I in the upper 
panel of the Figure. There may be an effect for higher 
metallicities, although that is where age differences be- 
come important as well. In the lower panel, we apply the 
age correction derived from Figure [5] and find that there 
is little if any residual metallicity dependence. Again, 
the principal driver of the trends appears to be age. 



Fig. 8. — Examination of structural effects on T e fm- We plot 
Y e FM using <TpM for the entire cluster sample vs. various struc- 
tural measures. In the upper panels, we plot T e fm uncorrected 
for the age relationship from Figure \5\ and in the lower panels 
we correct for the age relationship. Clusters older than 10 Gyr 
are represented by small open symbols and younger ones by filled 
symbols. The larger open circles again highlight those clusters 
with measured tr's, although we use ofm for all clusters. We plot 
Tjjtjk vs. a dimensionless measure of the central potential, Wo 
(leftmost panels), vs. concentration, c = log(r ti( ; a ;/r core ) (center 
panels), and vs. the distance of each cluster to the center of its 
pare nt galaxy (rightmost panels). That v alues of Wo and c come 
from McLaughlin & van dcr Marcl (2005) for their preferred Wil- 
son models. The magnitude of any correlations are significantly 
smaller than the age effects visible between the lower and upper 
panels. A small number of clusters (< 5) lie beyond the plotting 
limits, which are chosen to bring out details of the distribution of 
the bulk of the clusters. 
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Fig. 9. — The dependence of T e ,FM on cluster metallicity com- 
piled by McLaughlin & van dcr Marcl (2005) from a set of refer- 
ences. T e j?M is evaluated using (TpM f° r all clusters (small open 
circles). Clusters for which (jobs i s available are highlighted with 
large open circles. We recalculate T e fm using ctqbs and plot 
those results for comparison as filled circles. In the bottom panel 
we correct for the age dependence of T £ ^fm inferred from Figure[5] 
The age-related effects are significantly larger than any correlation 
with metallicity. 



3.5. Physical Sources of Apparent Scatter: Conclusions 

On the basis of the correlation between scatter and age, 
and the agreement between the inferred sense of the effect 
and the qualitative expectations from stellar population 
models, we conclude that age is the dominant physical 
source of scatter. However, our data do not allow us 
to eliminate all possible structural differences that might 
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Fig. 10. — Direct comparison of two ways to evaluate T e as a 
function of age. We plot the results using the stellar mass-to-light 
ratio from stellar population models, T t mo( j (small open circles), 
and from the FM using Equation (1) and <tfm> "^e,FM (filled small 
circles). Large open circles highlight clusters with <JQ B g. The ages 
for the values of T e p^j have been arbitrarily shifted by 0.2 for 
clarity. The two methods of estimating T e qualitatively track each 
other with age, but differ in the overall amplitude of the change 
from the youngest to oldest clusters. 

play a role, nor do they statistically demonstrate the su- 
periority of age as a driving parameter over the other 
parameters now available. The peculiarities of the sam- 
ple, only six clusters with measured a at age ~ 10 8 yrs 
and the remainder all at old ages, preclude application 
of statistical tests such as principal component analysis 
or partial correlation coefficients as a means to identify a 
dominant physical effect. However, this shortcoming can 
be remedied with reasonable observational effort. 

4. ESTIMATING T* 
4.1. Comparing Y^fm to Stellar Population Models 

The variation in T e , which for stellar clusters is equiv- 
alent to Y*, vs. age has a natural explanation in terms 
of stellar population evolution. In Figure [10] we com- 
pare estimates of T e based on the FM, T £t FM, to those 
ba sed on the model values of T», T t mn rf, tabulated 
by iMcLaughlin fc van der M arcl (2005) using their pre- 
ferred Bruzual-Charlot models with a Chabrier initial 
mass function (IMF). The sense of the behavior of T e 
with time is consistent for both T £i fm and T*, mo d sug- 
gesting that we are indeed seeing the effects of stellar 
evolution on T e . However, the quantitative agreement is 
poor, with T e j?M overestimating T* iTno d at young ages 
and underestimating it at old ages. 

The existence of a significant difference between 
T*. mo d and T Ct pM, particularly at old ages, is inde- 
pendent of the particular model or initial mass func- 
tion, among those currently in general use. We 
compare T e ,FM to the full range o f T* .m od pre- 
sented by Mc Laughlin fc van der Marell (|2005t) in Fig- 
ure 1111 We find similar results after comparing to 
20 07 Bruzual-Charlot models, which are not included 
in IMcLaughlin fc van der Marell (pOOST ). The magnitude 
of the differences between Y» imo£ 2 and T Ct pM is similar 
to that of the differences among the stellar population 
models themselves. It is therefore possible that there are 
problems with the stellar popul ation models at this level. 
iConrov. Gunn. fc White! (|2008|) explore the uncertainty 
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Fig. 11. — Comparison of T» jmo[ j and ~T £t FM for different stellar 
population models. Pane ls represent (a) Bruzu al-Charlot model s 
[IBruzual fc Chariot) 12001) with Chabrier IMF (ICharbrierl I2003T ): 
(b) Bru zual-Charlot models with Salp eter IMF; (c) PEGASE 
models HFioc fc R occa-Volmcrangc 1993) with Chabrier IMF; (d) 
PEGASE models with Salpeter IMF; (e) PEGASE models with the 
Chabrier globul ar cluster IMF E harbricr 2003); and (f) PEGASE 
models with the lKroupa. Tout, fc Gilmorel ^1993) IMF. The sense 
of the deviation for the oldest ages (highest values of T e ) is con- 
sistent among all varieties of stellar models. The stellar models 
consistently over predict T e relative to the FM estimates for these 
older clusters. 

in stellar population models, the bulk of which comes 
from modeling intermediate and low-mass stars near the 
end of their lives where they become far more lum inous 
but for short periods of time. iTonini et~a l. (2008) show 
that including the thermally-pulsating AGB phase can 
increase the emission from an intermediate age galaxy 
by 1 magnitude in the if-band. The level of the discrep- 
ancy suggested by our FM analysis is consistent with 
these potential problems in T* tmo d- Even at the young 
(low T e ) end of things, there are clear differences among 
the models, with those using the Salpeter IMF producing 
results more in line with those from the FM. Further dis- 
cussion of the details of stellar population models is be- 
yond the scope of this paper, but the FM provides a new 
constraint. The FM-derived values of T* suggest that all 
currently popular stellar population model variants over- 
pred ict the stellar mass -to-light ratio of old populations 
(see IMieske et all [2008, for an alternate approach that 
results in the same conclusion). 

Stellar population models that provide accurate stel- 
lar mass-to-light ratios are manifestly important for a 
variety of uses. However, star clusters on the FM pro- 
vide a new way to estimate T e for large sets of galax- 
ies. Currently, the FM zero point (C in Equation (1)) 
is calibrated using the results for T e from detailed dy- 
namical studies of sev eral tens of early-type galaxies 
(|Cappellari et all [2001 . If we knew, from stellar mod- 
els, the values of T„, and hence T e , for star clusters, we 
could use them instead to calibrate the FM. Once the 
FM is accurately calibrated, one can use it to solve for 
T e for any galaxy with measured V, r ei and I e . 

4.2. An Independent Prediction of T* 
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Fig. 12. — We equate T e j?k with Y* by assuming that star 
clusters are dark-matter free. In the upper panels we plot the 
inferred value of Y„, T« fm, vs. age and color. The curves are 
the by-eye fits to the data, with coefficients given in Tables [2] and 
\3\ The fit vs. color is valid for B — V < 0.75 and the data imply 
a constant value of Y», 1.13, for colors redder than 0.75. The 
bottom panels show the trends vs. metallicity. We can see what 
types of populations are not constrained by the current data (e.g., 
young, metal-poor populations) and for which the fitting functions 
fail (e.g., old, metal-rich clusters). 

We apply the FM to derive empirical relationships be- 
tween parameters such as main sequence turn-off age and 
color, and Y». Figure IT21 shows our crude fits to the age- 
T„ (Table [2]) and (B-V)-T* (Table E2 the fit is valid for 
B-V < 0.75, for B-V > 0.75 use T* = 1.13) relation- 
ships. Neither 1-D relationship appears satisfactory for 
the full range of clusters, and certain populations, such 
as young and metal-poor populations, are completely un- 
constrained. More sophisticated and well-constrained fits 
will be possible once more data become available for in- 
termediate age clusters, providing an alternative to T„ 
estimates that depend on stellar populations models. 

5. SUMMARY 
We find the following: 

1) On average, star clusters fall along the extrapolation 
of the Fundamental Manifold (FM) defined by spheroidal 
galaxies. Their larger scatter is consistent with the pro- 
portionally larger uncertainties in their velocity disper- 
sion measurements. Even extreme clusters such as Pal 
5 and 14 fall on the relationship when precise velocity 
dispersions are available. 

2) Individual clusters are offset from the FM in propor- 
tion to their age: older clusters fall on the manifold and 
younger clusters do not. Unfortunately, the sample of 
clusters with measured velocity dispersions does not in- 
clude many clusters with log (age [yr]) < 10, so we apply 
the FM relationship to infer a and T e simultaneously for 
the remainder. This procedure works faithfully for those 
with measured a. 

3) Aside from the trend with age, deviations from the 
FM do not correlate measurably with metallicity, dynam- 
ical state, or galactocentric distance. 

4) The sense of the evolution of mass-to-light ratio 
with age is as predicted by stellar population synthesis 
models, but currently popular mod els fail to reproduc e 
the estimates from the FM (see also|Micskc ct al. 2008). 
Specifically, the models over-predict the mass-to-light ra- 
tio of old, metal-rich populations. We provide empirical 



formulae with which to estimate the stellar mass-to-light 
ratios of certain cluster populations derived from the FM- 
estimated mass-to-light ratios. 

All stellar systems that we have ex amined in de- 
tail so far — brightest cluster galaxies (Zaritsky et al. 
I2006al). various galaxy types and lum inosity classes 
(jZaritskv et~aT1 l2008t IMieske et ail 120081) . and now star 
clusters — satisfy a simple, 3-parameter (2-dimensional, 
but not planar) scaling relation. For star clusters, de- 
viations from the scaling law are age dependent, which 
is qualitatively consistent with the expectation for pas- 
sive stellar evolution of the population. Quantitatively, 
however, this description fails. We suggest, but have not 
definitively demonstrated, that this disagreement arises 
from systematic errors in the stellar population models 
that lead to an overestimation of the V-band mass-to- 
light ratio of old stellar populations. The sense and 
magnitude of the effect are consistent with conclusions 
reached b y others investigat i ng stellar evolut ion models 
in detail ([Tonini et ail 120081: IMaraston et al.ll2006D and 
at the very least suggest that a factor of two uncertainty 
remains i n modeling stellar populations. As stressed by 
other s (cf. IMaraston et a l. 2006; Conrov. Gunn. fc White! 
2008), this uncertainty propagates in complicated ways, 
particularly for composite populations of unknown age 
distributions, in analyses of galaxy evolution. 

We conclude that there are no evident structural differ- 
ences, in either stellar spatial distribution or kinematics, 
between galaxies and star clusters beyond those of scale 
and those captured in the mass-to-light ratios within r e . 
Star clusters appear to be baryon-dominated versions of 
highly-compact, low-mass galaxies. As such, we are not 
surprised by the difficulties encountered by others in dis- 
tinguishing clusters from ultra-compact galaxies. If there 
are wholesale differences between the formation process 
and/or evolution of star clusters and galaxies, they do not 
manifest themselves directly in the gross physical prop- 
erties. We conclude that even if significant differences 
exist between the populations, due to the existence or 
lack of a dark matter halo, simple measurements will not 
reveal them. Projections of parameter space in which 
these populations are separate are deceptive and should 
not be taken to indicate that these objects are discrete 
families of gravitationally bound stellar systems. 

Although this lack of evident distinction is a set-back 
for attempts to differentiate the populations and probe 
potential differences in their formation and evolution, it 
is an advantage to those efforts to build a simple set of 
models for galaxies. Star clusters can provide an anchor 
for models of galaxies, empirically defining the stellar 
mass-to-light ratios or, alternatively, the zero point of 
the FM. The latter would enable us to use the FM to 
determine the mass-to-light ratio within r e , to a level of 
precision comparable to that obtained with detailed dy- 
namical modeling or strong gravitational lensing model- 
ing, for any galaxy with known V, r e , and I e . 

We return to the question of the defining character- 
istics of a galaxy. Among the stellar systems that we 
study, we have not identified a criterion that we consider 
a discriminating feature. Despite the popularity of the 
dark matter halo criterion, there exist "galaxies" that 
may not have a dark mat ter halo (tidal dwarf galaxies; 
iBarnes &: Hern quist 1992) or, at least, a d ominant dark 
matter halo (ultra-compact dwarf galaxies: IMieske et al.l 
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2008). In fact, from state-of-the-art optical data alone, it 
is difficult to demo nstrate that even giant ellipticals con- 
tain dark matter ([Cappella ri et al.| [2006). An alterna- 
tive, and more easily observable, classification is the rel- 
ative concentration of luminous baryons to dark matter. 
In such a scheme, systems composed solely of baryons 
are not distinct from those with dark halos and whose 
observed kinematics are dominated by baryons. While 
this scheme makes it harder to connect observations with 
dark matter simulations, the reality is that, for systems 
whose baryons are extremely concentrated, we cannot 
determine whether a halo exists or what its properties 
might be. 

Our primary empirical result is that t he existing de- 
scription of the Fundamental Manifold (Zaritsky et al. 
2008) is valid over an even larger parameter range than 
previously demonstrated and that now includes stellar 
clusters. The existence of a FM for the entire family 
of stellar systems implies 1) that the objects are virial- 
ized and 2) that there is limited variation in the overall 
structure. While the former is unsurprising, the latter 
implicitly places strong limits on how the systems form 
and how the luminous baryons settle. By extending the 
validity of the FM, or any scaling relation, over a larger 
range of parameter space and a wider set of objects, we 
are demonstrating that the physics that constrains the 
resulting baryon and dark matter distributions in stel- 
lar systems is more general than previously thought. In 
other words, whereas one might not have been surprised 
by the implication of the Fundamental Plane that all el- 
lipticals have related structural properties, it is certainly 
a puzzle to find that spirals, dwarf spheroidals, and now 
stellar clusters do too. 

While the structure of one stellar system is not a sim- 
ple scaling of another, there exists a transform function 
that relates the properties of stellar systems. Specifically, 
we have presented a fitting function for T e that is valid 
across all stellar systems. The generality of this trans- 



form function implies that 1) the structure of one type of 
stellar system is not a distinct result of a process partic- 
ular to that set of systems, such as violent relaxation for 
elliptical galaxies, but is instead the result of an inter- 
play of all processes responsible for the generic settling of 
baryons in gravitational potential wells, 2) the existence 
of the transform function for T e , and hence the nature of 
baryon settling, is independent of whether the system is 
embedded within a dark matter halo, and 3) differences 
in initial conditions or stochastic events during the life- 
times of stellar systems do not invalidate the transform 
function. By including stellar clusters, which we presume 
are devoid of dark matter, we demonstrate that the na- 
ture of the transform function for T e depends more on 
the baryon physics than the source of the gravitational 
potential. The presence or absence of dark matter im- 
prints no clear feature in the structure or kinematics of 
stellar clusters versus galaxies. By including stellar clus- 
ters, which we presume formed mostly in single, short- 
lived, dramatic fashion, we demonstrate that the nature 
of the transform function for T e is not altered radically 
by differences in formation history. The extension of the 
FM over a larger parameter space is in reality the ex- 
tension of the statement that physical processes rather 
than initial conditions or subsequent stochastic events 
dominate the final outcome in the formation of stellar 
systems. 
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TABLE 1 

FM-Predicted Cluster Velocity Dispersions 



NT 

1\ ill I1C 


&PRED 
[km sec ] 


&OBS 
[km sec ] 


Ref 


AMI 


0.6 






ARP2 


0.7 






PAL3 


1.2 






PAL4 


1.1 






PAL5 


0.7 


0.9±0.2 


Dehnen et al. (20041 


PAL12 


0.6 






PAL14 


0.5 


0.38 ± 0.12 


Jordi et al ("20091 



TABLE 2 

T*-(log(age)) Coefficients 



Order Number 1 2 3 4 5 6 



Coefficient 2282 -1736 545 -90.5 8.37 -0.409 -0.0082 



TABLE 3 

X, - (B - V) Coefficients 



Order Number 12 



Coefficient 



0.31 



-0.55 2.2 
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APPENDIX 

Appendix A 

Va r ious efforts have been made to establish a fun damental plane for globular clusters (Diorgovski 1995; McLaughlin 
2000; Barmbv et al.ll200"7t iPaquato &: Bertinll2008f ) . However, by necessity, this plane is not the same as that for ellip- 
ticals, whose Fundamental Plane (FP) requires T e oc <7 Q /f . T e for intermediate mass systems is already comparable 
to T* and cannot physically dec rease further for systems with comparable surface brightnesses but lower velocity 
dispersions, like globular clusters ()Korm cndv 1985|). 

A recent study d iscussed a unification of all "spheroidal" systems from globular clusters through elliptical galaxies 
([Forbes et alJl2008|). While prevailing generally, they found that dwarf spheroidals did not fit their general formalism 
and that the FM ([Zaritsky was not valid down to globular clusters. We disagree with the latter statement 

and so we investigate the nature of the discrepancy. 

We use their data, which they made publicly available with their paper, to remove any sample differences. Their 
claim regarding the FM is based on the nature of the empirical scaling relatio n for T e . In pa r ticula r, they note that 
the globular clusters do not follow the "U" -shaped pattern described for T e bv lZaritskv et al.l (|2008l ). 

Because their data include if -band magnitudes, our fitting formula for T e , which is calibrated to /-band values of 
T e , is not optimal (color differences become large). Instead we refit using their data and obtain a version of T e jn 
appropriate for the i-T-band. The central issue of the FM is whether the family of stellar systems can be described 
with a two-parameter relationship. The coefficients in the expansion of T e ja, or even the functional form, are not 
sacrosanct, and curr ently have no detailed physical justification. The result of the new fit is shown in Figure 1131 As 
iZaritskv et al.l (|2008f) stressed, the most salient feature of this Figure is not that the data straddle the 1:1 line, because 
this arises by construction (although it does imply that our functional form for T e ju is an acceptable one), but rather 
the low and uniform scatter across the full range of systems (i.e., clusters, dwarf galaxies, and gi ant galaxies). 

In the right panel of Figure[l3j we explore the crux of the argument against the FM presented bv lForbes et al.1 {2008). 
In particular, they noted that star clusters, with low values of T e , do not follow the upward trend in T e established 
by Local Group dSph galaxies. The difficulty with this argument is that it is based on a projection o f the d ata onto 
the (cr, T e ) axes. In the right-hand plot, we show the value of T e j it for each object in the lForbes et al.l (|2008|) sample. 
By construction, because we use the values on the manifold, there is zero scatter in this distribution, and yet in this 
projection there appears to be significant scatter. In fact, the clusters, which populate the lower left branch, appear 
to be a distinct population from the dSphs, which lie directly above the clusters and extend upward along the T e 
axis. Therefore, even points distributed on a surface with zero scatter can, in projection, appear to comprise different 
populations of objects because the surface is not uniformly populated. As such, one cannot conclude from plots such as 
that shown on the left of the Figure that there is no single scaling relationship among the populations. It is therefore 
critical to consider the behavior of Y e ju in terms of both a and I e . We find that the FM fits the characteristics of both 
high and low surface brightness objects with low (< 10 km sec -1 ) velocity dispersions. In particular, in our Figure Q] 
the FM was derived excluding the star clusters and UCDs, and yet both populations fall within the uncertainties of 
the extrapolated relationship. 

The apparent disagreement with iTollerud et al.l (|2010[) stems from their aims and approach. Their aim was to 
identify the parameter with the greatest influence in producing the full range of spheroid properties along the scaling 
relationship for spheroids. They correctly focused on velocity dispersion as the dominant parameter. They proceeded 
to derive a one-parameter scaling relationship that they could then tie to halo mass. Because they use only one 
parameter, they are unable to simultaneously fit the range of stellar systems with low a (i.e., dwarf spheroidals and 
star clusters/UCDs). Given their interest in tying stellar systems to dark matter halos, they naturally accepted 
the fit that was consistent with the dSphs, which are clearly embedded in dark matter halos, and argued that the 
clusters/UCDs are a different family distinguished by their lack of dark matter halos. That is certainly one viable 
interpretation given the state of the da ta for UCDs. However, we argue that 1) there is some evidence for UCDs 
having dark matter (Mies ke et al.l 120081 ) . 2) at the high-mass end, the UCDs overlap nicely with the low-mass end of 
the galaxy distribution (see Figure 1), and 3) disk galaxies also provide evidence for a two-parameter scaling relation. 
Our differences are therefore ones of interpretation regarding whether it is appropriate to consider all stellar systems 
together or to distinguish among them. Our preference for a single descriptive formalism for all stellar systems leads 
us in one direction, their preference for a 1:1 association of stellar systems with dark halos leads them in another. 
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Fig , 13. — The Forbes ct al. (2008) sample and the FM. On the left, we plot an edge-on projection of the FM using the Forb es et al.l 
(2008) sample of galaxies and star clusters. We refit for the function T e f it because they provide i-C-band, rather than optical, data. We 
have not attempted to account for age effects here. That the clusters show scatter comparable to the galaxies and are not significantly 
displaced from the relationship demonstrates that this sample of star clusters does satisfy the FM even prior to accounting for possible age 
effects. On the right, we place each of their objects on the fitting function surface to demonstrate how projections of the distribution of 
objects on this infinitesimally thin surface can create large apparent scatter. Judging whether a sample of objects satisfies the FM requires 
examining the "edge-on" projection, on the left, rather than the scatter between any two parameters alone. 



